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Abstract 

We consider the strong decays of the lowest-lying radially excited baryons, 
in SU(6) language the states comprising the first excited 56-plet. Assum- 
ing a single-quark decay approximation, and negligible configuration mixing, 
we make model-independent predictions for the partial decay widths to final 
states with a single meson. Masses of unobserved states are predicted using 
results from large- N c QCD, and the momentum dependence of the one-body 
decay amplitude is determined phenomenologically by fitting to observed de- 
cays, so that the baryon spatial wave functions are not assumed. We point 
out that comparison of these predictions to experiment may shed light on 
whether the Roper resonance can be consistently interpreted as a three-quark 
state. 
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I. INTRODUCTION 



In this note, we consider the strong decays of the lowest-lying radially excited baryons. 
These are the states comprising the first excited 56-plet, in SU(6) language. We make 
predictions for the decay widths of these states that are independent of any specific model 
for the binding potential or spatial wave functions. 

A reason for taking up studies of excited baryons now is that we anticipate new results 
from the CLAS Collaboration at the Thomas Jefferson National Accelerator Facility. Not 
all the states in the radially excited 56-plet are yet discovered, and many of the measured 
decay widths of the observed states have large uncertainties. A set of predictions for the 
decay widths of unobserved states will be of great use to experimenters. Of necessity, we 
also provide predictions for the masses of the unobserved states. To this end, we present a 
derivation of the Giirsey-Radicati mass formula [0 from large- N c QCD that we believe has 
not appeared explicitly in the literature. 

Another and more specific reason for our analysis is to shed light on the Roper resonance, 
the iV(1440). The Roper is the first excitation of the nucleon with the same J p = (l/2) + 
quantum numbers. The most direct explanation of the Roper is as a three-quark radially 
excited state. However, certain of its apparent peculiarities have led to plausible sugges- 
tions that it may be a hybrid state (a state whose lowest significant Fock component is 
three quarks plus gluonic excitations) or that it may be a cross section enhancement that 
does not correspond to any resonance [Q. An argument against interpreting the Roper as 
a three-quark state has been that its calculated mass is too high in quark models with 
one-gluon-exchange, i.e. spin-color, interactions Q. Recently, it has been found that the 
observed mass is perfectly consistent with the three-quark picture if mass splittings are due 
to spin-flavor interactions instead ||. Moreover, there is evidence of a Roper signal from a 
lattice calculation using a three-quark source operator, although uncertainties in the mass 
determination as yet preclude saying if it is lighter than its negative parity counterpart ||. 
In the present work, we assume the Roper resonance is a three-quark state embedded in the 
excited 56-plet; confirmation of the decay widths predicted here will therefore support the 
three-quark interpretation of this state. 

The value of a model-independent analysis should be clear: We are free of any as- 
sumptions regarding the binding potential or quark spatial wave functions. We do make 
assumptions about the dominance of single-quark operators and about configuration mixing 
that we now consider in turn. 

First, we assume that the decays of interest proceed via a single quark interaction vertex. 
We may write the decay operator in the same SU(6) tensor product space in which we define 
the baryon spin-flavor wave functions: 



Here G* a = S\T® is an SU(6) generator, where S\ and are the generators of SU(2) 
and SU(3), respectively, k l is the meson three-momentum, and ir a represents a meson field 
operator. The asterisk in Eq. ( |1 . 1|) indicates that that the operator acts on the spin-flavor 
wave function of the excited quark. Here we adopt a convenient picture of the baryon state 
that follows from large- N c QCD, namely, that the baryon consists of a 'core' of N c — 1 ground 
state quarks generating a collective potential in which a single quark is excited. Note that 
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the complete symmetry of the 56-plet spatial wave functions implies that matrix elements 



of the operator in Eq. (|1 . 1|) are proportional to matrix elements of the same spin-flavor 



operator summed over all the quarks in the baryon state, 

G m = £ S^T". (1.2) 

quarks a 

Thus, the operator in Eq. ( |1 . 1|) works equally well in parameterizing the decays of radially 
excited baryons for N c = 3, where the state may be described more realistically as a collective 
excitation of all three quarks. Matrix elements of this operator may be written 

TOIfteffl*^)) = f(k)tf(B f \GP a \Bi) (1.3) 

where k — \k\ and f(k) is a function that parameterizes the momentum dependence of the 
amplitude. This includes momentum dependence that originates from the underlying quark- 
level vertex, as well as from the overlap of the baryon spatial wave functions. The baryon 
states shown represent only the spin-flavor part of the wave function, and are nonrelativis- 
tically normalized. Use of the single-body decay approximation is strongly motivated by its 
success in describing decays of the orbitally excited, 70-plet baryons, both in the algebraic 
SU(6) approach |7J, as well as more recent large- N c effective field theory analyses pj-TTTIj. In 



the latter case, both one- and two-body operators have been shown to contribute to 70-plet 
decay matrix elements at leading order in the 1/N C expansion. Nonetheless, comparison to 
data for both strong decays |J and photoproduction amplitudes [fL(J reveals that the two- 
body operators are phenomenologically irrelevant. It is not unreasonable to assume that 
similar dynamics may play a role in the strong decays of other excited SU(6) multipletsQ. 

Our second assumption is that mixing between states in the excited 56-plet and other 
SU(6) multiplets, i.e. configuration mixing, can be neglected. In the case of the 70-plet, 
successful large- N c effective theory descriptions of the strong decays ||, photoproduction 
amplitudes |10 |, and the nonstrange mass spectrum M all neglect mixing with other SU(6) 



multiplets. In the present case, one might worry that the nearest observed states with ap- 
propriate quantum numbers for mixing, members of the positive-parity 70-plet, are only 
between 200 and 300 MeV heavier than the states of interest, so that a significant effect 
cannot be precluded. Fortunately, evidence from explicit quark models suggests that config- 
uration mixing involving the excited 56 is also small []12| |, and even smaller than one would 



expect from 1/N C arguments |TB[. We simply adopt this as a working assumption, without 
wedding ourselves to a particular quark model or choice of baryon spatial wave function. 



II. ANALYSIS OF OBSERVED STATES 

We first determine the functional form and normalization of f(k) by considering the 
decay modes of observed members of the excited 56-plet, the iV(1440), A(1600), A(1600) 
and the £(1660). We will refer to this multiplet as the 56' henceforth. (The RPP refers to 



Note that the suppression of two-body operators in the present case implies that F-wave decay 
amplitudes should be highly suppressed. 
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the 56' as the (56,0^).) Our input values are shown in Table |, and have been extracted 
from data in the Review of Particle Physics (RPP) fl4 |. The one standard deviation errors 
on the baryon masses were taken to be half of the corresponding mass ranges given in the 
RPP, and partial widths were computed using the quoted full widths, branching fractions, 
and associated errors. The values for f(k) corresponding to each observed decay mode are 
shown in Fig. |l|. Since we wish to model the function f(k) for momenta from ~ 100 to 
~ 500 MeV, we require only that we find a simple functional form that works well within 
this range. The most successful result is f{k) = (2.8 ± 0.2) /k, which has a \ 2 P er degree 
of freedom of 1.1. Note that the functional forms a + kb, a/k - 5 , a/k 1 - 5 lead to a x 2 P er 
degree of freedom of 1.9, 1.6, and 2.8, respectively. It is somewhat remarkable that such a 
simple functional form can account for the data in the momentum range of interest, though 
admittedly, the experimental uncertainties are large. It is not inconceivable, for example, 
that improvement in the data could lead to a preference for the less singular-looking linear 
fit, but given the present errors, the difference between a linear and the 1/k fit makes little 
difference in our decay predictions. Continuing in our spirit of model independence, we do 
not try to ascertain the origin of f(k), but explore how far we can go in predicting decay 
modes for unobserved states. 



III. MASS PREDICTIONS 

We use the Giirsey-Radicati (GR) formula to predict masses for the unobserved members 
of the 56-plet: 

M = A + BN S + C[I(I + 1) - N 2 /A] + DJ(J + 1). (3.1) 

Here J (I) is the baryon spin (isospin) and N s is the number of valence strange quarks. The 
GR formula predicts spin- and flavor- dependent mass splitting for completely symmetric 
SU(6) baryon multiplets in terms of the four parameters A, B, C and D; in our case, 
these parameters are determined from the masses given in Table |. While the GR formula 
was originally a conjecture QXj , and later 'derived' by assuming SU(6) symmetry and a 
specific set of mass operators transforming in symmetry-breaking representations of small 
dimensionality f[5" |, we begin by showing how Eq. ( |3.1| ) may be obtained more rigorously 



from large- A^ c QCD. This observation follows from the discussion in Ref. ||16|1 , but is not 
presented there explicitly. In a large- N c operator analysis for the 56-plet, we may write the 
general mass formula 

S 2 „ S l G i8 S 2 T S 2 T 8 T 8 

M = ail + a 2 — + ea 3 T + ea 4 — — h +ea 5 — — — h e a 6 



N c N c N 2 N c 

+ e 2 a 7 ^^ + e 3 a 8 ^^, (3.2) 

where S, T, and G are the spin, flavor, and spin-flavor generators of SU(6) (with a sum over 
quarks left implicit), and e ~ 1/3 ~ 1/N C parameterizes the size of SU(3) breaking. Notice 
that the eight order-one coefficients, ai . . . a 8 , completely span the space of observables, 
namely the eight baryon mass eigenvalues. Assuming that the baryon states have spin, 
isospin, and strangeness of order one in the large A^ c limit, we may discard all but the first 
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four terms if we choose to work only up to subleading order. Acting on a large-A c baryon 
state 

S 2 = S{S+1), (3.3) 

T 8 = (N C -3N S )/Vl2 } (3.4) 
and for totally-symmetric spin-flavor wave functions 

^G 18 = — ^ [3I{I + 1) - S{S + 1) - 3N S {N S + 2)/4] . (3.5) 

These identities imply that the only effects of the discarded terms up through order 1/N C 
are redefinitions of the first four coefficients, Thus we may write the 

mass eigenvalues 

a' 

M = a[N c + j^S(S + 1) + ea' 3 (N c - 3N S ) 



y C 



+ e^^- [31(1 + 1) - S(S + 1) - 3N S (N S + 2)/4] + (9(l/iV c 2 ) . (3.6) 



Eq. ( |3.1| ) then follows from a trivial redefinition of the coefficients. In modern language, the 
SU(6) breaking representations that lead to the GR mass formula are precisely those that 
arise at leading and subleading order in the 1/N C expansion. 

Using the observed masses in Table |, we may solve for the coefficients in Eq. ( |3.1|) . In 
units of MeV, we find 

A = 1406.3 ± 31.3, B = 195.0 ±43.0, 
C = 15.0 ±38.1, D = 43.3 ± 46.0 , 

from which we predict the following states: 

-'(1825 ±98), S*'(1790± 192), H*'(1955 ± 196), ft'(2120 ± 234) . 

IV. DECAY WIDTH PREDICTIONS 

From Eq. (|1.3|) , the decay width is given by 

where Mj and Mf are the masses of the initial and final baryons and 

|2| 2 = E \( B f\ G ja\Bi)\ 2 , (4.2) 

where the sum is over final state spins and isospins. The quantity \Q\ 2 may be obtained 
either by use of symbolic math code or by using 
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(Bf\G ja \Bi) = (Bf(If,af,Sf,Tnf)\G ja \Bi(Ii,a h Si,mi)) = Q 



h 
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(4.3) 



and evaluating the left-hand side for one particular final state. Here, a and m stand for 
isospin and spin projections, and the sum in equation (|4.2[ ) is over a/, a a , m,f, and j. 

We have already discussed fitting the function f(k) to the measured decays, and how to 
use the Giirsey-Radicati mass formula to obtain the masses of the unobserved 56' states. 
Using the results of the fit with f(k) oc 1/k, we predict the remainder of the strong decays 
of the 56'. The results are shown in Table |J, along with the masses we use for the initial 
baryons. 

Table [H] does not quote the uncertainties for the widths. Numerical uncertainties come 
from two sources. One is the uncertainty in the function f(k), and the other is the uncer- 
tainty in the masses. The uncertainty in the fitted f(k) leads to about ±15% uncertainties 
in the widths. Further width uncertainty induced by mass uncertainty can be large, as dis- 
cussed below. Of course, once the state is found, the widths in Table |I| can be recalculated 
easily and accurately with the correct mass. 



V. DISCUSSION 

The uncertainty in predicting the width of the 56' decays is dominated not by uncertainty 
in the matrix elements but by uncertainty in the masses of the states. Four of the 56' masses 
are measured and four are not. We must predict the unmeasured masses, and do so using the 
Giirsey-Radicati mass formula. The accuracy of masses predicted from the Giirsey-Radicati 
formula is ultimately limited by the the approximations inherent in its derivation, and is 
easily estimated in a large-iV c scheme. Currently, however, uncertainties in the predicted 
masses are dominated by uncertainties in the measured masses, and they are not small. 

Table || shows the predicted partial width of the various decay modes, using the central 
values of the predicted masses. To help see how the mass uncertainties affect us, consider 
the equal spacing rule that follows from the GR formula. The decuplet spacing is given by 

3 1 

Amity = 2 m ( A ') - 2 m< ^') ~ m(N') = (165 ± 110) MeV. (5.1) 

(For comparison, the equivalent prediction for the ground state decuplet is Am 10 = 139 MeV, 
to be compared to 153, 148, and 139 MeV for the A-E*,E*-H*, and mass splittings, 
respectively.) Then, as one example, the predicted central value for the Q' mass is 2.12 GeV, 
and this leads to a sizeable decay width for Q' — > S*K. However, the uncertainties allow 
the Q' to lie below the threshold for this decay. 



Moreover, a simple weighted averaging of the published results that the RPP JTJj uses as 
the basis of their mass estimates for the A' and E' gives somewhat different central masses 
and much tighter uncertainties than they conservatively quote. One finds 

m(A') = (1600 ± 7) MeV 

m(E') = (1672 ± 5) MeV. (5.2) 
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(The error limit on the A' mass includes a scale factor S = 1.4, in accordance with procedures 
described in the RPP narrative.) Using these masses together with the previous N' and A' 
masses yields 

Am w = (114 ± 20) MeV (5.3) 

This reduced spacing leads to a number of decays listed in Table [TT] being kinematically 
forbidden. 

Unfortunately, there is not much one can do theoretically to reduce the uncertainties 
originating from the 56' masses. Table H simply presents the most reliable predictions we 
can make with the current data. As unobserved states are discovered, and their masses 
measured, one can easily revise the widths in Table [II]. To reiterate, the relevant formula 
for our kf(k) = const. = a fit is 

kMf 9,^,9 , 

where M i; Mf, and k are defined already, the matrix element sum \Q\ 2 is given in Table ||, 
and a is 2.8 ± 0.20, based on the current data. 

The function f(k) is not taken from any model, but fit to the limited data with the 
fulfilled hope that there would be a good fit using a simple form. The form one would 
expect from an atomic physics calculation is rather different from the one we found. For 
the hydrogen atom, in a nonrelativistic calculation, f(k) would be proportional to a matrix 
element of exp(ik ■ r) between 2S and IS wave functions. Further in a hydrogen atom, 
the wavelength of the outgoing radiation is long compared to the Bohr radius, and the 
small k result that f(k) oc k 2 is valid. This leads to a very small width and the famous 
metastability of the hydrogen 2S state. For particle decay, the wavelength of the outgoing 
meson is comparable to and often smaller than the size scale of the baryon states. One 
should have no expectation that a low k approximation will work. This also applies to 
barrier factors in the non-met ast able case. 

A model or an eventual ab initio calculation for the baryon states will produce some 
definite f(k) that may or may not look analytically like our form, but to fit the data must 
not be numerically dissimilar to our result. 

To summarize, we have without reference to models, but using the limited available data, 
studied the strong decays of the lowest-lying radially excited baryons. We assumed, justified 
by previous studies, that one-body operators dominated and that configuration mixing could 
be neglected. Simple forms for the one momentum-dependent function give a good fit to 
existing data, and allow prediction of 22 additional decay modes for this multiplet. 
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FIG. 1. The function f(k). The best fit corresponds to / = 2.8/k, with x 2 /d.o.f = 1.1. 
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TABLES 





Mass (MeVJ 


Decay 


fc (JVleVJ 


Partial Width (JVleVj 


7V(1440) 


1450 ± 20 


Ntt 


397 


227.5 ± 67.3 








143 


OY.5 ± 30.5 


A(1600) 


1625 ± 75 


Ntt 


512 


61.25 ±31.5 






Att 


301 


192.5 ±76.0 


A(1600) 


1630 ± 70 


NK 


343 


33.75 ± 25.2 






Stt 


336 


52.5 ±51.3 


£(1660) 


1660 ± 30 


NK 


405 


24.0 ± 20.0 



TABLE I. Input values used in the determination of f(k). 
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Decay 


A //" /A A~A A 

iV2j(Mev ) 


k (MeVJ 


\r\2 


1 (Me V J 


A 7/ A T f 


1450 


397 


25/16 


164 


iV — ► A7T 1 


140U 


1 A Q 

14o 


Q 
I 


on 

y9 


A — > iV7T 1 


1625 


CIO 

512 


1 /I 

1/2 


60 


A / A „ t 
A — * A7T 1 


lozo 


3U1 


zo/lo 


145 


A — ► TV A 1 


1630 


343 


9/8 


91 


A — ► Z/7T ' 


1 con 


oof; 
OOD 


O //I 

3/4 


ID 


A — > 2j 7T 


1630 


188 


3/2 


97 


VV A7 TV' t 


1 r* n c\ 

1660 


405 


1/24 


4 


S' — ► A7T 


1 r* n c\ 

1660 


/i on 

439 


1/4 


on 

30 


Zj — ► Zj7T 


IbbU 


385 


2/<3 


f 


2j — > L 7T 


1 n c c\ 

1660 


O 1 O 

218 


1 /o 

1/3 


25 


c — ► AA 


1825 


403 


1/16 


6 


i — / \ -\ TV" 

C — > LA 


~i oo r 

1825 


323 


25/16 


1 O A 

134 


"/ ^-*" 


1 oo r 

1825 


A OA 

420 


1/16 


8 




1825 


240 


1/2 


A 1 

41 


2j —> Jy K 


1790 


519 


1/6 


1 O 

18 


L — > AA 


1790 


214 


5/6 


50 


Zj — ► A7T 


1 too 

1790 


for 

535 


1/4 


O A 

34 


2J — > S7T 


1790 


484 


1/6 


22 


L — > L 7T 


1790 


338 


5/6 


oo 

89 


L — ► Lr? 


1790 


196 


1/4 


13 


c — > AA 


i nrr 

1955 


C O C 

525 


1 / A 

1/4 


O 1 

31 


c — ► EA 


1955 


461 


1/4 


29 


S*' -> E*i? 


1955 


239 


5/4 


86 




1955 


520 


1/4 


36 




1955 


358 


5/16 


36 




1955 


269 


1/4 


19 


n' -► HA" 


2120 


506 


1 


128 


0' -► E*K 


2120 


274 


5/4 


101 



TABLE II. Decays of the radially excited 56' into ground state 56 plus meson. Seven of these 
decays are measured (indicated by a dagger), and the other 22 widths are predictions. If the mass 
of the initial state is known, the uncertainty in the width from our fitting procedure is ±15%. 
Uncertainties in the mass induce further and sometimes large uncertainties in the width. These of 
course are greatly reducible once the mass of the state is measured, as discussed in the text. 
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